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We present results of micromagnetic simulations of thin ferromagnetic rings undergoing magne- 
tization reversal. This geometry is one of few examples in micromagnetics in which the transition 
states have been found analytically in a ID model. According to this model, at low fields and large 
ring sizes, the energetically preferred transition state is a localized magnetization fluctuation (in- 
stanton saddle). At high fields and small ring size, the preferred saddle state is a uniformly rotated 
magnetization (constant saddle). In the first part of this paper, we use numerical micromagnetic 
simulations to test these predictions of the ID analytical model for more realistic situations, includ- 
ing a variety of ring radii, annular widths and magnetic fields. The predicted activation energies 
for magnetization reversal are found to be in close agreement with numerical results, even for rings 
with a large annular width where the ID approximation would be expected to break down. We find 
that this approximation breaks down only when the ring's annular width exceeds its radius. In the 
second part, we present new met ast able states found in the large radius limit and discuss how they 
provide a more complete understanding of the energy landscape of magnetic nanorings. 



I. INTRODUCTION 

The magnetic properties of thin ferromagnetic annuli 
have attracted attention due to their potential applica- 
tions in magnetic random access memory: the absence in 
such geometries of edges or corners to nucleate magneti- 
zation reversal leads to greater stability against reversal 
than in other simply connected thin film elements. 

There are several ways in which ferromagnetic an- 
nuli may be used as memory elements, differing in the 
(meta) stable magnetization configurations that represent 
a single bit. In all of these the magnetization lies com- 
pletely within the plane and its configuration is smooth 
everywhere. One such pair of configurations is of oppo- 
site chirality, i.e., clockwise or counterclockwise circula- 
tion of the magnetizationiii^i^; another is the so-called 
"onion states^i^i^i^i^i^, where there is a net total magne- 
tization along a direction in the ring's plane. In the for- 
mer (latter) case, a circumferentially (uniformly) directed 
magnetic field can be used to set the magnetization con- 
figuration. The minimum energy configurations depend 
on the strength and direction of the magnetic field as well 
as the relative dimensions of the ring with respect to the 
exchange length. 

Few analytic solutions have been obtained for the rate 
of thermally induced reversal in micromagnetic problems, 
even in relatively simple geometries. For thin annuli 
under the influence of a circumferential magnetic field, 
however, the lowest energy transition (or saddle) states 
have been found analytically^^ in a one-dimensional ap- 
proximation and studied numerically in the full three- 
dimensional problem^. 

The Kramers theory of reaction rate^ can be used 
to compute the typical lifetime of a given state when 
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Figure 1: Ring geometry showing the coordinates. The cur- 
rent runs along the axis of the ring out of the page. 



there are several minimum energy configurations. The 
general form for the rate F of thermally induced tran- 
sitions between two minima in the limit of low noise 
is given by the well-known Arrhenius formulaic F = 
Fo exp{—AE/kBT), where the prefactor Fq is usually in- 
dependent of the noise strength and depends only on the 
shape of the energy landscape close to the extremal states 
relevant to the transition. The activation energy, AE, 
equals the energy difference between the transition state 
and the metastable state, thereby determining the sta- 
bility of the latter. This is an important figure of merit 
for memory devices. 

In this paper, numerical micromagnetics are used to 
test the predictions of the analytical theory of Martens et 
al^i^ for thermally induced transitions between states of 
opposite chirality in a ID approximation to the annular 
ring. The simulations were made for a variety of mean 
radii, annular widths and magnetic fields. 

The geometry under study and accompanying relevant 
parameters are represented in Fig. [TJ The magnetic ma- 
terial is in the shape of an annulus of mean radius R, an- 
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nular width AR and (in the third dimension) thickness t. 
A current / running along the axis of the ring produces 
a circumferential external field H(r) = {I/27rr)0. The 
ring is composed of a soft isotropic ferromagnet (e.g., 
permalloy) with saturation magnetization Mg and ex- 
change length A. In all of the geometries considered, 
the aspect ratio /c = t/i? <C 1, giving rise to magne- 
tostatic forces that constrain the magnetization to lie 
in the plane of the ring (M^ = 0)^^. A magnetization 
configuration can therefore be completely described by 
0(6>,r), the angle the magnetization at a given radius 
makes with the unit vector lying along the tangent to 
the circle with that radius: M(l9,r) = {M^^My^M^) = 
M,(sin((/)-6>),cos((/)-6/),0). 

II. MODEL 

Our starting point is the Landau-Lifshitz- 
Gilbert (LLC) equationiiii^ 

^ = -|7|M X Heff - X (M X Heff), (1) 

where 7 is the gyromagnetic ratio and a is the (phe- 
nomenological) damping constant. The effective mag- 
netic field, Heff = — Vm£^, contains all (external and in- 
ternal) fields and is the variational derivative of the total 
micromagnetic energy 

E[M{r)] = M^J dV|VMp 

(2) 

+ f /r3 d^r\\/U\^ - /io fa ^'^He • M. 

The three terms above correspond respectively to the ex- 
change energy, demagnetization (or magnetostatic) en- 
ergy, and Zeeman energy, with the (small) magneto- 
crystalline anisotropy term neglected. (The last of these 
can be easily included, but for the materials and ge- 
ometries considered here, it is typically overwhelmed by 
the much larger shape anisotropy arising from the de- 
magnetization term.) Here Q is the volume of the ring, 
A = y/2A/ (/io^s ) is the exchange length (where Ms, 
the magnitude of the magnetization, is assumed to be 
the same everywhere, and A is the exchange constant). 
He is the applied external magnetic field, and |VMp = 
|VM^P + |VM^P + |VM^p. The "magnetostatic poten- 
tial" U, arising from long-range dipole-dipole interactions 
within the magnetic material, satisfies V^t/ = V-M (and 
suitable boundary conditions in the interfaces between 
media), which can be derived through Maxwell's equa- 
tions. Our simulations involve numerical integration of 
the above set of equations. 

The extremal states of a quasi- II) ferromagnetic ring 
(i.e., AR <C so that the external magnetic field does 
not vary significantly with distance from the center of 
the ring) in a circumferential magnetic field have been 
analytically obtained^^. The solutions found there ap- 
ply in the thin-ring limit: k = t/R <^ 1 and (A/i?)^ ~ 



{t/R) \ ln{t/R)\. Under these conditions the second term 
on the RHS of ([2]) separates into three main terms (and 
a number of smaller ones) : a term which extracts a large 
energy cost when the magnetization does not lie com- 
pletely within the plane of the annulus; a local surface 
term (the shape anisotropy, which favors alignment of the 
magnetization with the tangential direction at the inner 
and outer ring radius); and a nonlocal bulk contribution. 
Analysis of these terms finds that the bulk term is small 
compared to the surface term and can therefore be ne- 
glected^^. In Sec. [lV]we test these conclusions for more 
realistic geometries by computing numerically the total 
demagnetization energy and comparing it to the (analyt- 
ically computable) local surface (i.e., shape anisotropy) 
term. 

In the ID approximation the total energy reduces to 

(3) 

xA2 /^"[(^If )2 + sin^ - 2/icos(/.]d0, 

where the parameters £ and h are the scaled circumfer- 
ence and field: 

^=fV'2-(^)|l-a)l 

(4) 

and He is defined below. The first term in the integrand 
is the exchange energy, the second the shape anisotropy 
(i.e., the surface term arising from the demagnetization 
energy), and the last is the Zeeman term. 

Given an external magnetic field that is circumferen- 
tial and points everywhere in the counterclockwise di- 
rection, there are two states that are local minima of 
the energy: a stable magnetization configuration (ground 
state), which is everywhere aligned with the external 
field, and a met ast able state that is everywhere antipar- 
allel to the field (i.e., circumferential and pointing ev- 
erywhere in the clockwise direction). He corresponds in 
eq. [4] to the magnetic field at which the metastable con- 
figuration becomes unstable. 

There are also two relevant unstable stationary config- 
urations (i.e. saddle states). These are defined by the 
angle that the magnetization direction makes with the 
circumferential direction at each point in the annulus; 
i.e., (/) is a function of the angle (as shown in Fig. [T]), 
(pi^hifi). In the limit of low noise, reversal of the magneti- 
zation occurs through the lower energy saddle state. One 
of these corresponds to a global rotation of the magneti- 
zation in which (j) is independent of 9] we therefore label 
it the "constant saddle", and is denoted (})h • The constant 
saddle favors the exchange and Zeeman energies at the 
expense of the demagnetization energy. The second sad- 
dle state is a localized fiuctuation of the magnetization 
and we therefore refer to it as the "instanton saddle", and 
denote it by (I) h/{0) {denoted as instanton saddle ini^). 
This state favors the demagnetization energy at the ex- 
pense of the exchange and Zeeman energies. 
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Which of these two saddles is energeticahy favored de- 
pends on the apphed field and the ring siz e. When the 
scaled field h is smaller than — {27i/lY the instanton 
saddle has a lower energy than the constant saddle; oth- 
erwise, the constant saddle is lower in energy. Fig. 2 ofl^ 
shows the phase boundary between the two activation 
regimes as a function of h and I. 

We simulated the dynamics using the analytical so- 
lutions as our initial configurations for the two saddle 
states: 



(t)h cos~^{-h) 

for the constant saddle, and 

1 , , /OK(m) 
Dh£ = 2 cot"^ ( 1? dn I ^—\m 



(5) 



(6) 



for the instanton saddle. Here dn(-|m) is the Jacobi ellip- 
tic function with < m < 1, and K(m) is the complete 
elliptic integral of the first kind^^. The parameter m sat- 
isfies. 



m 



2K(m) 



(7) 



and i9 is defined by. 





2 — mh — Y^^^^^^^^^^^"^^^0^^-^^ 



2m — 2 + mh + w?h? + 4(1 — m) 



(8) 



For sufficiently small rings [i < 27r) , the instanton sad- 
dle does not exist (essentially, the variation of the mag- 
netization, which is of order the exchange length, cannot 
"fit" onto the ring). In this limit m ^ and the in- 
stanton solution reduces to the constant saddle. As the 
ring becomes larger the parameter m increases from 
to 1 monotonically with i. For £ ^ 27r the m becomes 
numerically indistinguishable from 1. In this limit the 
instanton saddle configuration is given by: 



2 tan 




^cosh(-Vr 



(9) 



We classify rings according to their i values as small 
{£ < 27r), medium {£ > 27r), and large {£ > 27r). Phys- 
ically, i characterizes the ratio of the magnetostatic to 
exchange energies. A medium ring {i ^ 27r) has a scaled 
circumference close in size to a domain wall in the mate- 
rial. The saddle configurations of each regime are shown 
in Fig. El 

The activation energies within the ID analytical model 
can be calculated using [Sjcfji^). For the constant 
saddle they can be analytically computed: AE = 
^oM^t^R\\n{t/R)\{l - h)^ . 



Figure 2: Saddle configurations computed analytically from 
the ID model of Martens et al for small, medium and large 
ring sizes, a) For small rings the magnetization reversal goes 
through the constant saddle state; b) when i > 27v the tran- 
sition is via the instanton saddle; c) as the relative size of the 
ring increases the fluctuation in the instanton saddle occupies 
a smaller fraction of the ring. 



III. METHOD 

We studied thin nanorings by running simulations on 
the model of Sect. [Til using the publicly available pack- 
ages OOMFF and Nma g^^i^^ . These packages effectively 
simulate the dynamics specified by ([T]) and ([2]) at zero 
temperature; i.e., all runs start from an initial configura- 
tion and run downhill in energy. 

Our initial states were the instanton and constant sad- 
dles described by (O and ([6]), which provided starting 
points that were guessed to be relatively close to the ac- 
tual saddles. The system subsequently relaxed to the 
actual saddle states, which turned out to be remarkably 
close to analytical solutions. We describe below how this 
was determined. 

Depending on the starting state, the system will evolve 
to one or the other (meta) stable state, i.e., either the 
clockwise or the counterclockwise magnetization configu- 
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ration. In order to find the actual 2D saddle numerically 
(recall that the magnetization is forced by the magne- 
tostatic energy term to lie in the plane of the ring), we 
introduce a new field value, denoted /it, an "effective field" 
for which the state (j)h,i behaves as a saddle state. In de- 
termining ht^ two criteria are used. First, (l)h^i must be as 
nearly a stationary state as the numerics allow, i.e., the 
initial time derivative of the total micromagnetic energy 
should be close to zero (lim/^^^/^^ "WL^o ^ 0~). Sec- 
ond, the state (l)h^i should mark the boundary between the 
basins of attraction of each (meta) stable state (i.e., for 
h < ht the system evolves to a clockclockwise state, while 
for h > ht, it evolves to the counterclockwise state). In 
the following sections, we show that these criteria are sat- 
isfied in a variety of rings with different exchange lengths 
and annular widths. We also show how the model even- 
tually breaks down when the width of the ring becomes 
very large. 

The procedure can be summarized is as follows. For a 
given initial (j)h,i we find the appropriate ht by a bracket 
and bisection iterative process. We set the initial config- 
uration fix the external magnetic field at the value 
he and allow the system to relax. If the final state is the 
metastable (clockwise) configuration we increase he by 
an amount She] if the final state is the stable (counter- 
clockwise) configuration we decrease he by She • We then 
start a new run and reduce She by a factor of 2. As She 
decreases, the total relaxation time increases due to the 
slow dynamics at the start of the simulation, providing 
evidence that the initial configurations are approaching 
the true saddle states. We iterate until we reach a nu- 
merical uncertainty of Sh = 6 x 10~^. 

The Nmag simulations were run using a mesh consist- 
ing of 20963 volume elements, 15154 surface elements, 
and 7594 points with an edge length of average 3.89 nm 
and standard deviation of 0.7. (The quality distribution 
of the mesh was 1.66% below 0.6; 9.11% between 0.6 and 
0.7, 57.28% between 0.7 and 0.8, 32.05% between 0.8 and 
0.9, and 0.01% above 0.9.) The cell sizes for the simula- 
tions in OOMMF were selected so that they were smaller 
than the exchange length for each of the regimes studied. 
Different values of i for a given geometry were studied 
by changing the exchange constant A (cf. below ([2])) and 
keeping the ring dimensions constant. This changes A 
and therefore ^, the mesh can be adjusted to speed up 
the simulations. 



A. String Method in Rare Events 

The String Method^^ is a recently introduced numer- 
ical procedure for calculating transition energies and 
paths within the context of large fiuctuations and rare 
events. It is useful to find the path connecting two sta- 
ble configurations Ma and Mb, through a curve ^ with 
minimum energy. The obtained path corresponds to the 
reversal trajectory in the limit where the precession term 
of ^ is negligible compared to the damping term. This 



curve ^ satisfies 

Vm^^(0 = Vm^(0 - [Vm^(0 • t\i = (10) 

where t is the unit tangent of the curve ^. The curve ^ is 
found by guessing a parametrized path ^(0) = {M(a) G 
[0,1], M(0) = Ma,M(1) = Mb} and evolving it in 
"time" according to: 

dtM{a) = -VME^{M{a)) + A£ (11) 

The second term is added to enforce a particular 
parametrization; it does not alter the actual evolution 
of the curve. It is convenient to rewrite this equation as: 

dtM{a) = Hi-f^{M{a)) + Xt (12) 

For numerical purposes the path ^ is discretized with 
N-\-l points betwen MAand Mb- After each iteration of 
(p!2|) with an Euler forward algorithm the magnetization 
vectors are renormalized to Mg. In Sect. [VTlwe use this 
method to find the barrier between two states connected 
through a transition state. 



IV. MEDIUM SIZE NARROW RING 

We consider a narrow ring of medium reduced circum- 
ference £ (i.e., the parameter m not close to 1) with 
X/R < 1, t/R <C 1 and (X/R)^ - {t/R)\\n{t/R)\ 
{A = 3.2 X 10-^^ J/m, AR = 40 nm, R = 200 nm, 
t = 2 nm, Ms = 8 X 10^ A/m). With these values, £ and 
He are 12 and 73.9 mT, respectively. We first test numer- 
ically whether the surface, or shape anisotropy, term is 
in fact the main contributor to the magnetostatic energy, 
as required for the vailidity of the analytic solutions to 
hold^°. Using OOMMF and Nmag the total demagneti- 
zation energy were obtained for different values of h and 
compared to the values of the surface term (second term 
in the integrand of (|3])). The results of this comparison 
are presented in Fig. [31 which 

shows that the numerical computation of the total de- 
magnetization energy gives nearly the same dependence 
on h as the contribution from the shape anisotropy alone. 
This provides numerical support for the approximations 
used to arrive at the ID analytical solutions of^^ and 
confirms that the bulk magnetostatic term, neglected in 
the analytic model, is indeed not important. 

We also compare the total demagnetization ener- 
gies computed using OOMMF for different cell sizes in 
Fig. [3t). Three suitable methods to calculate the demag 
energy are available in the OOMMF package; ConstMag, 
3dSlab, 3dCharge. Constmag calculates the average de- 
magnetization field assuming the magnetization is con- 
stant in each cell; 3dSlab uses a demag field obtained 
from blocks of constant charge; and 3dCharge uses con- 
stant magnetization to calculate the in-plane component 
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Figure 3: Total demagnetization energies, using different cal- 
culation methods, compared to analytical results for the shape 
anistropy term alone {A = 3.2 x 10~^° J/m, R = 200 nm , 
t = 2 nm and AR = 40 nm). a) Comparison of OOMMF 
using the ConstMag method at cell size 1 nm and Nmag with 
theoretical predictions vs. h. b) Demagnetization energy vs. 
cell size using three different methods for computing this en- 
ergy in OOMMF, all for h = 0.2. The soUd line is the analyt- 
ical computation of the shape anisotropy term alone. 



of the magnetic field, and constant charges to calculate 
the out of plane demagnetizing field. As seen in the fig- 
ure, the consistency between different methods of calcu- 
lation improves as cell size is reduced; and the numerical 
results approach that of the ID model. 

Once ht is obtained following the method described in 
Sect. lIIIt the saddle state is numerically obtained and the 
activation energy is thereby determined from the differ- 
ence between each saddle state and the metastable state. 
Fig. [4] displays curves at ht = 0.21 for each of the two 
saddle states used as initial configurations. 

As can be seen in the figure, after a very short tran- 
sient the system arrives at a configuration in which the 
energy stays almost constant for an extended period; this 
indicates that the initial ID analytical solution is close in 
both energy and its geometrical configuration to the true 
2D saddle. Eventually, the saddle state decays into one 
of the two stable configurations. The activation energy is 
easily computed this way (as seen in Fig. [4]), and a glance 
at the figure confirms that for the applied field ht = 0.21 



0.2 0.4 0.6 0.8 1 1.2 

t(ns) 

Figure 4: Evolution of the total micromagnetic energy with 
time. Arrows show activation energy of each saddle. The 
configurations shown bifurcates at the fields given and have 
an initial slow evolution. The two criteria used to determine 
a saddle state. 



the instanton configuration has a lower activation energy 
than the constant saddle, as predicted theoretically^^. 

The method described in Sect. [TTll was repeated to 
obtain the behavior of AE as a function of ht. Us- 
ing this approach one can calculate ht ~ ^^"^ and 

^ E[^,,ht]+^E[^,,h:] _ _ ^^^^ 

two saddle configurations. Here, he] represents the 

numerical energy of a configuration under an applied 
field of magnitude he. The results are summarized in 
Fig. [5] together with the analytical predictions. 

Figs. [5^1 and[5t) show AE{ht) for the instanton saddle 
and constant saddle, respectively. From bottom to top 
the curves represent calculations in OOMMF for two dif- 
ferent cell sizes, an Nmag calculation and the analytical 
prediction. 

The Nmag results are closer to the analytical predic- 
tions than the OOMMF calculations. This is due to the 
fact that the curvature is more faithfully represented by a 
mesh of tethrahedrons in Nmag, whereas OOMMF rep- 
resents the ring with a square grid. The edges of the 
tethrahedrons can be arranged to follow closely the cur- 
vature of the ring, whereas the square cells edges will 
in most cases make a finite angle with the ring tangent. 
This results in a much larger error contribution to the de- 
magnetization energy in square grids in OOMMF than in 
Nmag. 

Fig. [St presents the Nmag simulation results for the 
activation energy. As predicted in Ref.^^, the instan- 
ton saddle configuration is preferred at lower fields. For 
a fixed i the activa tion energy curves are predicted to 
cross at hc{i) = ^/T^^J27^/ly . For higher fields, the 
constant configuration is the sole saddle state (this is in 
contrast to the low-field side, where both solutions ex- 
ist but the constant saddle has higher energy). Numer- 
ically, the field at which the crossover between the sad- 
dles occurs is somewhat lower than that predicted. This 
discrepancy arises because the theoretical model applies 
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Figure 6: Evolution of the total energy with time for the 
instanton saddle, with energies measured with respect to that 
of the stable state. The decay to the metastable state is shown 
as proceeding to the left of t = 0, and to the stable state to 
the right of t = 0. The ring dimensions are AR = 40 nm, 
i = 60. When he < ht the instanton saddle (II close to 
t=0) decays quickly into the metastable configuration (I at 
t=0.7 ns); domain walls are annihilated. When he > ht the 
two domain walls of the instanton saddle (III. close to t=0) 
separate (IV, t=1.5ns), move to the opposite side of the ring 
(V, t=5.4ns) and annihilate (VI, t=5.7ns). 



Figure 5: Comparison of the activation energies for the (a) 
instanton and (b) constant saddles. Data points represent 
steps of 0.1 in h. As the numerics are refined both by using 
Nmag and by reducing the cell size in OOMMF simulations 
the results approach the theoretical predictions, (c) Compar- 
ison for Nmag for the middle size, narrow ring {£=12, Ai?=40 
nm) of the corresponding curves shown in (a) and (b). 



to a strictly ID ring, whereas the simulations run us- 
ing higher-dimensional rings (2D in OOMMF and 3D in 
Nmag) . 



V. LARGE SIZE NARROW RING 

We now investigate an annular film with the same 
dimensions but different A (e.g., permalloy with A = 
1.3 X 10~^^ J/m); the rest of the dimensions are kept 
the same as above. Such a ring belongs to the large ring 
regime {£ = 60, m 1). As before, we can obtain a qual- 
itative understanding of the reversal process by following 
the time evolution of the micromagnetic energy; this is 
presented in Figs. l6|7[ Using the fact that an instanton is 
a superposition of two domain walls with opposite chiral- 
ity (Fig.fTOb)'^. the features observed in Figs. [6l [71 can be 
explained and an intuitive idea of the reversal mechanism 
developed. 

We begin by discussing the evolution of energy with 
time as shown in Fig. [6^. For < ht the instanton de- 
cays into the metastable state with a rapid decrease in 
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Figure 7: Evolution of the total energy with time for the con- 
stant saddle, with energies measured with respect to that of 
the stable state. The decay to the metastable state is shown 
as proceeding to the left of t = 0, and to the stable state to the 
right of t = O.The ring dimensions are AR = 40 nm,^ = 60. 
When he < ht the constant saddle (II at t=0.15ns) starts to 
decay into the metastable configuration (I at t=0.28 ns), it is 
possible to see that this decay is not uniform in all the ring, 
and some regions decay faster than others. When he > ht the 
constant saddle quickly develops domains (III at t=0.51ns). 
Counterclockwise domains expand at the expense of clockwise 
domains (IV, t=0.83ns), with annihilitation that produce a 
sudden decrease in energy (V, t=1.03ns) leaving a single do- 
main wall pair which continues moving (VI, t=2.52ns) until 
eventual annihilation; the slope during the last stage is very 
close in magnitude to the slope of Fig. [6] during domain wall 
propagation, the discrepancy is due to the difference of ap- 
plied fields. 



energy (seen to the left of t = 0). This corresponds to the 
annihilation of two transverse domains walls. The field 
in this case is not sufficient to separate the two trans- 
verse walls. When he > ht (right side) the energy of 
the system decays linearly with time towards the stable 
state and then sharply decreases. The linear time de- 
cay corresponds to movement of the transverse domain 
walls in opposite directions (due to the external field) at 
approximately constant speed as can be seen in Fig. [6l 
The final, sharp decrease in energy results from the col- 
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Figure 8: Total energies of the two saddle states for dimen- 
sions AR = 40 nm,-^ = 60 with respect to the clockwise con- 
figuration. 



lision and annihilation of the walls. Note that the mag- 
netostatic and exchange energies remain almost constant 
while the walls are propagating (roughly) independently 
of each other. The annihilation energies of the domain 
walls (indicated by the arrows in Fig. [6fci) are roughly the 
same magnitude. 

The slope of the E{t) curve during the propagation 
phase provides a measure of how fast the reversal pro- 
ceeds. For the particular damping parameter used, this 
can be used to estimate the wall speed from the last term 
of ® to be t; = )2^(f )where n is the num- 

ber of domain wall pairs present in the ring. In deriving 
this expression cj) = 0, tt for each domain and domain wall 
widths are assumed to remain constant during propaga- 
tion. A comparison to medium-sized rings shows that at 
fixed h the reversal time increases as A decreases for two 
reasons: the domain wall width is comparatively smaller 
and the effective scaled circumference i is comparatively 
larger (cf. Figs. [Hand [6]). 

The time dependence of the energy for the constant 
saddle case. Fig. [3 is seen to proceed in several steps 
where the energy decrease is gradual, punctuated by large 
changes in the slope dE/dt. These features can be ex- 
plained as follows. Fig. [8] shows that the activation en- 
ergy of the constant saddle at 0.17 is several times 
larger than the activation energy of the instanton saddle 
at the same field. With such a large activation energy, 
it is relatively easy to create several domain wall pairs 
along different (randomly placed) parts of the ring, each 
of them with an energy cost roughly equal to the activa- 
tion energy of the instanton. In the simulations, this pro- 
cess is modified by the discretization of the ring, and in 
an experimental setup it is possible that impurities and 
edge roughness might play a similar role. The abrupt 
changes in slope are due, as before, to the annihilation 
of domain wall pairs and hint at a richer spectrum of 
states that are stationary in the energy when the scaled 
circumference is large compared to the exchange length. 
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Figure 9: Topological defects of the XY model in the bulk 
(a,b) and in the presence of an edge (c,d), with winding num- 
bers of 1,-1, +1/2,-1/2 respectively. 



A. Multiple Wall Pair Spectra 

We have also gone beyond the work of Martens et al 
by finding numerically new stable states consisting of 
multiple domain wall pairs. When present, these states 
infiuence the time evolution of the system in a manner 
similar to that just described. In this section we discuss 
these new states and present a model to incorporate their 
effects on magnetization reversal. 

Fig. \TOti shows one of these locally stable double pair 
states. They can be described as a combination of topo- 
logical defects, in particular edge defects and domain 
walls, which we will now describe (cf. Fig. [9]). Bulk 
topological defects are vortex and antivortex singularities 
of the magnetization configurations with a net contribu- 
tion to the exchange energy^^. They are characterized by 
their winding number which is conserved over any con- 
tinous transformation of the magnetization. Close to an 
edge of the material the singularities become half-vortices 
with winding number ±1/2^^. Since i ^ oo for the limit 
m ^ 1, we can consider any small segment of the ring 
as a strip: the outer edge of the ring maps into the lower 
part of the strip. 

A transverse domain wall can be described as a com- 
posite of two edge defects of opposite sign^ at opposite 
sides of a ferromagnetic strip (see Fig. [TOk). In the ring, 
it is convenient to label such a domain wall using the 
sign of the topological defect on the inner side of the 
ring. The topological defects experience a "Coulomb-like" 
attraction or repulsion. Walls where the magnetization 
points in opposite directions (equal signs for same side 
edge defects as in Fig. [TOb ) experience repulsive inter- 
actions. The origin of repulsion arises from the mag- 
netostatic and exchange energies in the region between 
walls. Walls which are parallel to each other (with oppo- 
site signs in the same side edge defects, as in Fig. [TOb) 
experience attractive interactions. 

Any small fiuctuation of the magnetization initially 
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Figure 10: a) Transverse domain wall composed of two edge 
defects of opposite sign, b) Magnetic configuration equiva- 
lent to the instanton state in a stripe, c) A trapped domain 
between two antiparallel walls, d) Low energy metastable 
state configuration. The topological defect sign at the in- 
ner boundary determines if the domain wall pairs are stable 
with respect to the external magnetic field. With respect to 
the ground state (counterclockwise), this state has a winding 
number of zero, as do the instanton saddle, constant saddle 
and clockwise configurations. 



parallel to the ring (strip) edges would be a precursor to 
a double wall of this last type (Fig. [TOb) as illustrated by 
the profile of such configurations. The radial component 
of the magnetization in the ring (transversal component 
in the strip) edges could have any sign: the fiuctuation 
will have the same energy whether the magnetization tilts 
towards the inside or the outside of the ring (up or down 
in a strip). A reversal can be produced by a fiuctuation 
with equal probability for any of these two orientations. 

The domain between two walls will expand under the 
infiuence of a parallel external magnetic field (producing 
a repulsive pressure on the walls) , and contract under an 
antiparallel magnetic field (producing an attractive pres- 
sure on the walls). The balance between the interdefect 
interaction and the field determines whether a configu- 
ration is in stable or unstable equilibrium. For example, 
an instanton saddle configuration is equivalent to two do- 
main walls with opposite signs on their innermost defects 
which enclose a domain parallel to the field (cf. signs in 
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Figure 11: Two possible evolutions of a fluctuation with two 
instantons in different parts of the ring. (A) The magnetiza- 
tion relaxes to the metastable state only if the radial magne- 
tization components are parallel. (B) Otherwise, two trapped 
domains (360 degree walls) appear. 



Fig. [To] b). The field pushes the domain walls away from 
each other while their mutual interaction tends to bring 
them together. These opposing tendencies produce the 
unstable equilibrium which makes this configuration a 
saddle state. 

The opposite situation, in which edge defects repel and 
the enclosed domain is antiparallel to the field, produces 
a metastable state. The enclosed magnetic domain does 
not vanish because the half vortices experience an ef- 
fective repulsion: it is energically costly to unwind the 
topological defects. As a result this configuration is sta- 
ble, with an energy intermediate between the clockwise 
and counterclockwise configurations. 

We can now explain how the metastable state evolves 
into the state represented in Fig. [TOl i. If the ring is orig- 
inally magnetized clockwise, the energy density per ring 
segment of instanton fluctuation size is higher than the 
instanton energy. If the ring is sufficiently large, two in- 
stantons (one domain wall pair each) are produced (Fig. 
[TTj) . If both instantons point in the same radial direction 
the system evolves into the counterclockwise state (Fig. 
[TT]A). However, if the instanton fluctuations point in op- 
posite radial directions (Fig. HUB) the system evolves to 
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Figure 12: Spectra of large ring regime (m ^ 1). Lower 
metastable states are separated by the energy of a trapped 
domain configuration. Beyond the clockwise configuration the 
instanton saddle and the constant saddle are shown. The 
metastable configuration becomes unstable beyond /le = 1. 



the state shown with 4 domain walls . 

The domain wall pairs need not be at opposite sides of 
the ring for this configuration to be stable, as confirmed 
by displacing one of the wall pairs by several angles and 
waiting for the system to relax. 

It is interesting to summarize the new possible configu- 
rations and their total energies in the large ring case. The 
lowest state is the stable configuration, and there exists 
a series of metastable states separated from each other 
by the energy of a double wall trapped domain (Fig. [12]) 

With the exception of the single wall pair, all config- 
urations shown have a total winding number difference 
from the counterclockwise state equal to zero for a path 
that completely encloses the central hole of the ring. The 
single wall pair configuration has this winding number 
difference equal to one. While reversing the field will 
make the double wall system decay into the stable state; 
the single wall pair configuration cannot decay into the 
counterclockwise configuration. A trapped domain con- 
figuration using a single pair of domain walls has already 
been proposed for an MRAM device^. 



VI. WIDE RINGS 

Having verified numerically the conclusions and pre- 
dictions of the analytical model of Martens et al.^^ for 
narrow rings, we now proceed to test the limits of its ap- 
plicability with increasing annular width. Given the ID 
nature of the analytical mode, we expect that at some 
width the model should break down and its conclusions 
no longer apply. Surprisingly, this breakdown finally oc- 
curs at a larger width than initially expected. 

Increasing the ring width introduces two new effects 
that cause the analytical model to break down. First, 
it allows the external field to vary in magnitude appre- 
ciably as one moves along a radial direction. Second, it 
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Figure 13: Activation energy spectra for R = 200 nm, 
and a) AR = 40 nm, b) 100 nm, c) 200 nm. 



=12 



increases the relative magnitude of the (previously ne- 
glected) nonlocal bulk term of the magnetostatic energy 
with respect to that of the local surface term. 

In our simulations, AR is set to the values 100 nm 
{He 29.5 mT), 200 nm {He 14.7 mT) and 380 nm 
{He =7.8 mT). The central hole of the annulus is a few 
exchange lengths in diameter in the latter case. For this 
reason vortex-like "singularities" remain in the gap. 

Fig. [13] summarizes the activation energies for the 
widths considered. It should be noticed that as AR in- 
creases for fixed A, £ decreases and the annulus shifts 
away from the large ring (m ~ 1) approximation. We 
consequently discuss only middle size rings where i ^ 12. 
Fig.[T3]a,b,c shows that the main predictions of the model 
hold even for very wide rings: the configurations (j)h are 
saddles for certain ht, and for fields below he the in- 
stanton saddle configuration is preferred to the constant 
saddle configuration. 

We have found that the annular width must be in- 
creased to the extreme wide-ring limit AR ^ 2R in 
order for the model to fail. Its breakdown can be ob- 
served in the E{t) curve of Fig. [T4k . In this regime, 
there are still fields ht for which the dynamics bifurcates 
to either the stable or metastable state on either side 
of htj but at such fields it is clear from Fig. [T4k that 
{liiRh^ht ^It-o 7^ This indicates that the initial 
configurations chosen from the ID analytical solutions 
are no longer close to the true saddles. Instead of a long 
initial period of little change, we find instead relaxation 
to a state in which the central region of the ring is mag- 
netized circumferentially and the outer edge of the ring 
retains some memory of the starting configuration. This 
appears to be a new type of saddle configuration, which 
we call the relaxed state, and is shown in Fig. [T4l 

For simulations starting from the instanton saddle with 
parameters h = 0.1, ...,0.4 the system evolves to the re- 
laxed state. At higher values of h {h = 0.5, ...,0.8 for the 
instanton, and h = 0.7, 1.0 for constant saddle config- 
urations) the relaxed state does not satisfy the station- 
arity condition, but the bifurcation condition can still be 
satisfied with a particular ht. The absence of a plateau 
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Figure 14: a) Time evolution of the total energy, for both 
initial saddle configurations. At a given ht the dynamics do 
not slow arbitrarily as t ^ 0. b) Spectra of the approximate 
activation energies (measured close to the bifurcation in time, 
not at t = 0). c) The new metastable state of Fig. [15] is found 
at lower energies than either saddle. The ID saddle configu- 
rations (instanton center-up, constant center-down) and their 
corresponding relaxed states at t=0.8 ns for {ht right, h~ 
left); inside the ring the strong field aligns the magnetization 
clockwise . 



in the E vs. t curves makes the definition of the activa- 
tion energy somewhat more problematic, but it can still 
be defined by using as the energy of the saddle state the 
point at which the curves E[(j)h^ h~] and E[(j)h^ h^] sepa- 
rate. In this way approximate activation energy spectra 
can be determined, and the results are shown in Fig.[T4b. 

For values of h that correspond to the constant saddle 
(0.2,. . .,0.6) another feature of the breakdown of the ID 
model can be seen. The system relaxes to neither of the 
stable states considered so far (see Fig.[T5|). It evolves to 
a radially dependent state with counterclockwise orienta- 
tion at the inner edge of the ring and clockwise orienta- 
tion at the outer edge. This configuration is stabilized by 
the large inhomogeneity in the magnitude of the magnetic 
field as one moves outward along a radial direction. It is 
important to note that the energy of this state is lower 
than the energies of either initial configuration used (the 
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ID analytical saddle solutions), but is higher than either 
the clockwise or counterclockwise state. 

There exists a low barrier that prevents this configura- 
tion from relaxing into either of the counterclockwise con- 
figurations in the limiting fields of Fig. [HI We obtain an 
estimate of this energy by using the String Method with- 
out reparametrization^^ as described in Sect. Hill We 
start with sequences of 50 equally spaced configurations 
that connect each of the circumferencial configurations 
to this newly found metastable state along an straight 
line. The result of the relaxations are shown in Fig. [161 
Fig. Uni (a) shows the total energy of the string points 
after relaxation; the inset shows a very shallow energy 
barrier that prevents decay to the two lowest stable con- 
figurations. It is interesting to observe that this barrier 
almost disappears on the left at he = 0.22500 and on the 
right at he = 0.25625 explaining why this state is stable 
only for a very narrow band of field values. 

The physical origin of the two barriers is clear if the 
three components are studied separately (Fig. [161 b, c,d). 
For all graphs the origin corresponds to the state repre- 
sented in Fig. [m and the clockwise oriented magnetiza- 
tion is 50 steps to the left of zero, the counterclockwise 
configuration is located 50 steps to the right of zero. 

The Zeeman energy (Fig. [T6li) prevents the system 
from moving towards the clockwise configuration for all 
states along the string. It is maximum for the clock- 
wise configuration and decreases monotonically along the 
path. Any point of the trajectory is pushed to the right 
of the graph by the external field. 

Fig. [T6b shows the demagnetization energy with a 
sharp barrier that prevents the magnetization from point- 
ing perpendicular to the ring edges. At this barrier, the 
magnetization at the surface points radially outward (as 
in Fig. [T4h^ for the constant saddle). This produces 
a sharp barrier at this step of the path. At this point, 
the magnetostatic energy is the only energy term acting 
against the reversal of the magnetization from clockwise 
to counterclockwise orientations; other terms favors the 
reversal. The net effect of the demagnetization energy 
barrier is to favor configurations away from this barrier. 
Fig. [T6b represents the exchange energy. The exchange 
energy is minimal in the two circumferencial states. 

When these three interactions are considered together 
the stability of this state is understood: at lower fields 
the exchange and magnetostatic energy are balanced by 
the Zeeman energy. At large fields Zeeman and exchange 
favor a magnetization out of the ring's surface, when the 
shape anisotropy barrier is overcome, the system reverses 
suddenly into the counterclockwise configuration. 

The analytical model presumed the field to be con- 
stant in the radial direction as in the narrow ring case. 
Although the field inhomogeneity eventually renders this 
assumption invalid, the ID model is surprisingly robust 
and breaks down only in the extreme case just consid- 
ered. 



j/f f t f f / f 

Tt 1 1 t t t / / 



^ 4^ ^ IH^^^^^^^P ^ >r 

/ / I V \ \ \ 

^ i i \ \ 

^ / / A 



Figure 15: Magnetization configuration of new found 
metastable state for extremely wide rings. The stability of 
this state is a consequence of the highly inhomogeneous field 
in the center. 
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Figure 16: Energies of the configurations along the string 
after 100 iterations at two different applied fields. The cen- 
ter of the figure represents the magnetization configuration 
of Fig. [151 The clockwise configuration is at -50, and the 
counterclockwise is located at +50. Energies are a) total b) 
demagnetization c) exchange d) Zeeman . The inset on (a) 
shows the energy on a finer scale, illustrating the existence of 
an energy barrier. 



VII. CONCLUSION 

The ID analytical model of Martens et al}^ has been 
tested and confirmed using numerical simulations for a 
variety of ring sizes and external magnetic fields fields 
that more closely approximate realistic laboratory situ- 
ations. Although it was initially expected that the an- 
alytical model would apply only to narrow annuli, our 
simulations show that it is surprisingly robust, eventu- 
ally breaking down only in the extreme two-dimensional 
limit. 

By studying a large portion of the relevant parameter 
space, (A, Ai?), we have also found new saddle and stable 
states. These findings enrich our understanding of the 
energy landscape of ferromagnetic rings. 

Two limits present particularly interesting features: 
the large- narrow ring (Ai? <C R) and the extremely 
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wide ring. The large narrow ring allows for the appear- 
ance of multiple instantons at energies below the constant 
saddle configuration; their relative orientations and po- 
sitions determine the final magnetization configuration. 
We provide a topological analysis of these new configu- 
rations and predict a hierarchy of such states differing 
by the number of domain wall pairs in each. The in- 
teraction between these domain walls can be character- 
ized and understood through "edge defects" that compose 
them. Moreover, by following the evolution of the down- 
hill energy run from one of these states to the stable 
configuration, one can infer the propagation of domain 
wall pairs in the ring: sudden changes in energy indicate 
annihilation of domain wall pairs, while linear decrease 
of energy occurs during domain wall propagation. 

The ID model predicts extremely well the activation 
energy for magnetization reversal even for wide rings. 
Eventually, though, in the extremely wide regime limit, 
the 2D nature of the magnetic field becomes important 
and the ID approximation breaks down. In this regime, 
a new stable state arises in which the magnetization is 



radially dependent but independent of the angle. The 
String Methodic has been used to verify the existence of 
a barrier between this state and other stable states with 
lower energy, and is used to clarify the physical origin of 
this barrier, in terms of the various contributions to the 
energy from exchange, magnetostatic, and external field 
sources. 
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